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CMBR distortion concerned with recombination 
of the primordial hydrogen plasma 

E.E. Kholupenko 1 , A.V. Ivanchik 2 , and D.A. Varshalovich 3 

Ioffe Physical Technical Institute St. -Petersburg 

CMBR distortion concerned with recombination of the primordial plasma is calculated in frequency band from 
1 GHz to 100 GHz in the frame of the standard cosmological model for different values of cosmological density 
parameters: nonrelativistic matter density f2 ra and baryonic matter density fit, . Comparison of these results with 
observational data which will be obtained from planned experiments may be used for independent determination of 
the cosmological parameters Q m and Qt . 



1. Introduction 



The Universe was existing at the thermodynamical 
equilibrium during 100 thousand years after the elec- 
tron - positron annihilation. When the temperature 
of the Universe became low enough ( < 2 • 10 4 K, it 
is corresponding to redshift z < 7000) because of the 
cosmological expansion, the recombination of the main 
elements began (Fig. 1). 
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Figure 1: The recombination of the main elements: the 
solid line corresponds to the recombination of HII, dashed 
- Helll, dashed-dot - Hell 



This process leds to the birth of non-equilibrium 
photons. The most of them had survived until the 
modern epoch and have appeared as lines of the CMBR 
spectrum (Fig. 2). 
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Figure 2: Top picture: AJ is the brightness 
[erg/cm 2 /sec/Hz/ster] of the non-equilibrium radiation 
concerned with the primordial hydrogen recombination. 
Bottom picture: J [erg/cm 2 /sec/Hz/ster] is the CMBR 
brightness as a function of frequency. 



The measurement of the intensity and the form of 
these lines would give the information on the values of 
key cosmological parameters (e.g. f2b,f2 m ). 

The calculation of the distortion of the CMBR 
spectrum requires: 

1. The calculation of the level populations of the pri- 
mordial hydrogen atoms in the presence of a strong 
radiation (photon per mode r\ 3> 1 for the transitions 
investigated) which has the Planck spectrum at tem- 
perature decreasing with the cosmological time. 
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2. The modelling of the radiation transfer in the ho- 
mogeneous expanding environment. 

2. Physical model 

We used the standard cosmological model (see Tab.l). 
The temperature, the concentration of protons, and 
the Hubble constant as functions of redshift are: 



T = T (l + z), Ntot = N p0 (l + z) 3 



(1) 



H = iW^A + fim(l + Zf + VL re l(l + z) 4 (2) 



Table 1 . The standard cosmological model parameters 



Value description 


Symbol 


Value 


total matter 


Qtot 


1 


space curvature 


n k 





non-relativistic matter 




- 0.3 


baryonic matter 


n b 


- 0.04 


relativistic matter 




- 10~ 4 


vacuum-like energy 


n A 


- 0.7 


Hubble constant 


H 


70 km/s/Mpc 


radiation temperature 


To 


2.725 ±0.002 K 


proton concentration 


N p0 


- 2 • 10~ 7 cm" 3 



The characteristic values of the cosmological parame- 
ters in the recombination epoch are presented in Tab. 2. 
The knowledge of these values allows us to estimate 
the relative CMBR temperature distortion concerned 
with the non-equilibrium radiation at i — > k transition 
(hvik <C kgT) at frequency v ~ Vik/{^ + z) by the 
formula: 



AT 



1 



hvi, 



tonA Az k B T{z) 



N tot {z)0 ik 



(3) 



where Az is the duration of recombination epoch, 
9ik is an element of the matrix, so-called "matrix of 
efficiency of radiative transitions" (ERT-matrix was 
introduced by Bernstein et al. [1]). Their values are 
10~ 2 — 10~ 4 for i, k = 10 — 40 correspondingly (see [2]). 

Table 2 



Redshift 
Temperature 
Proton concentration 



z 
T 
N„ 



800 - 1600 
2200 - 4100 K 
100 - 700 cm" 3 



For example, according to formula the relative 
distortion of the CMBR temperature concerned with 
17 — > 16 transition (redshifted frequency ~ 1.14 GHz, 
9 ~ 5- 10~ 3 ) is AT/T = 2.5 • 10~ 8 . 

Formula (JSJ does not take into account the follow- 
ing facts: 

1. There are overlaps of the line profiles because the 
cosmological expansion leads to the formation of the 



wide line profile with the relative frequency width 
about 0.3. 

2. The variation of the recombination rate leads to the 
formation of a non-rectangular line profile. 

3. The elements of the ERT-matrix Oik depend on red- 
shift. They vary during recombination epoch within 
20%. 

The points mentioned above made us to calculate 
the formation of the CMBR distortion numerically. 

3. Initial equations 

3.1. Description of population behavior 

We considered the hydrogen atoms in the excited states 
with principal quantum number n > 2 ([1], [2]). The 
transitions into the ground state are considered sepa- 
rately because the optical depth for these transitions 
is much more than one ([3]). 

For the description of the behavior of the hydrogen 
atoms we used the quasi-stationary approximation of 
the kinetic equations averaged over the angular mo- 
mentum: 



oo 

E 

fc=l 



AkiXk - ^2 AikXi + onN e x v - PiXi - Ji = (4) 



k=l 



where Xi — Ni/N to t is the relative population of the 
state with principal quantum number n = i + 1 (sub- 
script i = l corresponds to the first excited state) , N{ 
is the population of state i (the number of atoms in 
state i per volume unit), N e is the free electron con- 
centration, x p is the free proton fraction, A ik is the 
coefficient of i — > k transition (An =0), 014 is the 
coefficient of recombination to state i , is the coeffi- 
cient of ionization from state i . Value Ji is the rate of 
uncompensated transitions from state i to the ground 
state of hydrogen atom. Value Ji is an independent 
parameter and to be described below. 

In this work, the free electron concentration N e (t) 
and the free proton fraction x p (t) are also independent 
parameters. These functions are the solution of the 
cosmological recombination problem. For the first time 
this problem has been solved by Peebles [3] , Zel'dovich 
et al. [4] for hydrogen primordial plasma. The problem 
of the recombination of the hydrogen-helium primor- 
dial plasma has been solved by Seager et al. [5] who 
used the model based on [3]. 

We have reproduced the results of investigations [3] 
and [5] (see Fig. 3). 

To solve equations Q the following symbols are 
convenient: X is the vector of the values Xi, B is the 
vector of the values hi = aiN e x p , J is the vector of 
the values Ji, Q is the transition matrix which has the 
elements given by the following expressions: 



A k i, i ^ k; 



\fc=i 



■0* 



(5) 
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Figure 3: The ionization fraction of the hydrogen plasma 
as a function of redshift 



These definitions allow us to rewrite equations (jlj : 

Q(t)X - J[t) + B{t) = (6) 

The solution of this equation can be presented in the 
form X = X° + AX , where X° is the equilibrium val- 
ues of the populations given by the Boltzmann distri- 
bution, AX is the corrections concerned with the un- 
compensated transitions into the ground state. Value 
X° is the solution of the equation 

Q(t)X° + Bit) = 0, (7) 

then AX is the solution of the equation 

Q{t)AX = J(t). (8) 

The explicit form of AX value is 

AX = Q- 1 (t)J(t) (9) 

The method mentioned above (equations0]-|nj) requires 
preliminary determination of J value. The exact cal- 
culation of this value demands of the joint modelling 
of: 

1. The recombination process. 

2. The population behavior, taking two-quantum de- 
cay of high-excited states into account. 

3. The Lyman-radiation transport for each line n — ► 1. 
In this paper we used the following approximation: 

all uncompensated transitions occur only from the first 
excited state due to two-quantum decay of state 2S and 
L a -quantum out of a line profile ( [3] , [4] ) . Thus we can 
write Ji — \ip\5\i {5i k is the Kroneker symbol). 

This definition allows us to write the following ex- 
pression for the population corrections: 

Axi = {Q~ l )a\x p \ (10) 



System (JSJ has been solved by the Gauss method with 
the choice of the leading element. The equilibrium val- 
ues have been determined by the Boltzmann distribu- 
tion. 

To show a convergence of calculations with a vari- 
ation of the level number, we used three models of 
the hydrogen atom: 40-level, 80-level, and 160-level. 
The convergence is quite satisfactory for our aim (see 
Fig. 7). The standard calculations are performed for 
80-level model. 

3.2. Description of radiation behavior 

In order to describe the radiation behavior, we used 
the kinetic equation for radiation: 

9v _ Hv Hm = 

dt dv 

2 oo i— 1 

= 8^2^*°* Y.J2( A * kX * ~ A kt x k )^ lk (v) (11) 

i=l k=Q 

where r]{t, v) is the number of photon per mode as 
a function of time and frequency, ^i k (v) is the line 
profile normalized by the expression J ^i k (v)dv = 1 . 

The function ^i k {v) has the maximum at fre- 
quency v — v ik and differs noticeably from zero in 
the frequency range [vd, v u ] , where v& ~ v ik — ZAvr) , 
v u ~ + 3Az^j , Avry — 3 • lO^'Vifc is the thermal 
line width. 

The solution of the equation Ijllfl can be presented 

as: 

oo i 

Tj(t,v) =Vo(t^)+J212 A ^(t^) (12) 

i=2 fe=l 

where rj (t, v) is the solution of the homogeneous equa- 
tion corresponding to equation i|ll|) , Ar\i k (i, v) is the 
partial solution of the equation: 

dArj lk _ Hi/ dA?jik _ 
dt dv 

c 3 

= - — ~N t ot{Ai k Xi - A ki x k )^ ik (v) (13) 

OTTV 

The property of the equation lllfl allows us to consider 
the equation for each line separately then to sum up 
the partial solutions. 

The problem of the behavior of the non-equilibrium 
radiation can be solved in two stage considered below. 

3.2.1. Formation of non-equilibrium radiation 

To research the formation of non-equilibrium radiation, 
we considered the frequency range [vd, v u ] close to cen- 
tral frequency of the line. The quasi-stationary 
approximation in this range ([3]) is: 

dAn s c 3 

-Hv— -ia = - — ^N tot (A ik x l -A kl x k )^ lk (v) (14) 
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where Ari? k (t,v) is the partial solution in the range 
considered. This solution depends on time paramet- 
rically. The value Ar)? k (t, v) at the lower limit of 
the range considered sets the boundary condition for 
the further evolution of the non-equilibrium radiation: 
^Vikit) = ^Vik( v d,t) ■ This value can be obtained by 
the formula based on the general solution of the radi- 
ation transport equation in the case of the interaction 
between radiation and matter at different tempera- 
tures ([6]): 

A Vik = (Vrn - f7o(fik))(l - CX P (-T lfc )) (15) 

where value t] m — l/(xki 2 /xik 2 — 1) describes the ra- 
diation of matter, t/q is the photon number per mode 
in the flow of the incident quanta. In our case rjo is 
given by the Planck distribution. The value Tik is the 
optical depth in the line i k ([7], [8]): 



„3 As -2 

c at ^ik i 1 „ 
Tik = z T^tot—rry-r^Xk 



(16) 



where A\ k is the Einstein coefficient of the spontaneous 
transition ([9], [10]). The formula l|l(jfl is obtained by 
means of the absorbtion coefficient formula of line i — > 
k ([11]). 

We can linearize expression <|15ll over the optical 
depth, if it is much less than one (t^ <C 1): 



and using l|lfi(l it can be written: 



(17) 



8nv; 



-N tot (z)H 1 (z){A ik x i -A kl x k ) (18) 



The following definition is convenient: 

Oik = Ji 1 (A ik x i - A ki x k ) 
It can be rewritten (in our model): 

Oik = \x p \^ 1 (A ik Ax i - A ki Ax k ) 



(19) 



(20) 



Taking into account the relation dz/dt = —H(z)(l+z) 
and expression (|20|l . formula (|18|) can be rewritten in 
the form: 



A 4 



87TJ/V 



-N tot (z)(i + z)-^e lk 

dz 



(21) 



This formula allows us to obtain the following (in the 
case hi>ik <C fc^T): 



AT 



(1 + ^ 



dx„ hu- 



lk 



dz ksT(z 



■N tot {z)0 ik (22) 



For the first time this formula has been obtained in the 
paper [1] in somewhat other form. 

We can estimate dx p /dz on its averaged value 
dxp/dz = 1/ Az. Then we obtain formula ©. 



The element of the ERT-matrix 0^ means the av- 
eraged number of photons emitted in line i — > k per 
one act of uncompensated recombination. 

Taking definition H20|l and expression p0(l into ac- 
count, the elements of the ERT-matrix are given by 
the formula: 



ik = (AikiQ-^n - AuiQ- 1 )^) 



(23) 



It shows that the elements of the ERT-matrix can be 
calculated independently with the solution of the cos- 
mological recombination problem. 

In paper [1] the calculation of the ERT-matrix was 
necessary because of the approach accepted. In paper 
[2] the ERT- values were the main result. To compare 
our result with that of paper [2] , ERT-matrix has been 
calculated in this work (Fig. 4, Fig. 5). The Burgin's 
results [2] have been reproduced with relative accuracy 
within 1%. 
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Figure 4: The ERT-elements as functions of the principal 
quantum number n of lower level: the solid line corre- 
sponds to transition with m = n + 1 , dashed - m = n + 2 , 
dashed-dot - m = n + 3 



3.2.2. Radiation transport in homogeneous 
expanding environment 

To investigate the radiation transport in the homoge- 
neous expanding environment, we considered the fre- 
quency range [0, v d ] . 

Due to the properties of function ^ik{v) , equation 
(|11|) in this range is: 



0_ 

di 



Arjik 



d 

Hv—A^ik 
av 







(24) 



where Arjik (i> v) is the partial solution with the bound- 
ary condition Ari lk {t,Vik) = Ar]f k (t) (taking v ik ~ v d 
into account). 
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Figure 5: The relative change of ERT-elements during the 
recombination epoch: from top to bottom the transitions 
are showed: 3^2, 11 -»• 10, 21 -> 20, 31 -» 30 



Changing from time- functions to redshift-functions, 
equation i|24|l takes the form: 

d d 
{l + z)— Ar) ik +v— Ar) ik = (25) 
oz ov 

with the boundary condition Arjik(z, Uik) = Ai]f k (z) . 
The solution of the equation (|25|l is: 

Arnk(z, v) = Ar)f k ((l + z)u ik /v - 1) (26) 

The number of photons per mode at frequency v 
in the epoch with redshift z is given by the formula: 

oo i 

V( z ,v) = Vo(z,v) + ^2^2 ^Vtk(z,f) (27) 

i=2 k=l 

The relative distortion AT/T = (T ex - T)/T of the 
CMBR temperature is calculated by the formula: 

Te*(z, u) = - h " (28) 

k B m(l/ri(z,v) + 1) 

4. Features of observation 

Many lines cannot be registered because the CMBR 
is contaminated by quasar radiation and IR-source ra- 
diation ([12]). For example, the strongest expected 
line L„ is invisible because of dust radiation. Taking 
into account this background, it is reasonable to cal- 
culate the recombination radiation concerned with the 
transitions with the principal quantum number over 
ten. The spectrum of this recombination radiation is 
located in the frequency range observed with the de- 
tectors for CMBR measurements. 



5. Results 

The main result of this paper is the relative CMBR 
temperature distortion as a function of frequency. 
This result is shown in Fig. 6. We calculated the 
CMBR distortion with the various sets of the cosmo- 
logical parameters il m and $!{,: {57 m = 0.3; h 2 ilb = 
0.014, 0.018, 0.022} and {tt m = 0.1, 0.3, 0.5; h 2 fl b = 
0.018} 




V, GHz 

Figure 6: The relative CMBR temperature distortion as 
a function of frequency. The curves correspond to h 2 flb = 
0.014; 0.018; 0.022 from bottom to top. All curves were 
calculated at Q m = 0.3 



The result of Jl m -variation is not presented graphi- 
cally because of weak dependence of CMBR distortion 
on £l m -value. This fact can be explained for the opti- 
cally thin lines (in our case r < 10~ 4 ) by formula (|21|l . 
It displays that, the CMBR distortion Aij depends on 
VL m by only the ionization fraction x p {z) which vary 
with Q m weakly. 

The fif, -variation shows that CMBR temperature 
distortion is directly proportional to the baryonic mat- 
ter density fit, (Fig. 6). 

The maximum value of the relative CMBR distor- 
tion in the observed range is less than 3 • 10" 7 . This 
value is to be observed at the low frequencies. 

To resolve the distortion corresponding to differ- 
ent values fif,, it is necessary to measure the CMBR 
temperature with absolute accuracy ~ 1 nK. 
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Figure 7: The distortion at variation of the level num- 
ber of hydrogen atom. The curves correspond to n max — 
40, 80, 160 from bottom to top. All curves were calculated 
at n m = 0.3,h 2 n b = 0.018 
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